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Abstract 

We introduce a new numerical technique to compute mass spectra, 
based on difference method and on a new gauge fixing procedure. We 
show that the method is very effective by test runs on a SU{2) lattice 
gauge theory. 
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Even if computations of the gluonic mass spectrum have reached recently 
very high levels of precision and reliability [1] looking for more effective phys- 
ical procedures to measure correlation functions is a very important task. 
The measurement of connected correlation functions from Monte Carlo pro- 
cedures is a very time consuming task, and the interesting signal comes from 
large separation, where the functions we want to determine are exponentially 
small. 

In this note we will continue in exploring the possibility of using con- 
tinuous updating schemes (Langevin-likc, as opposed to the discrete step 
Metropolis scheme) and to exploit such continuity to define connected cor- 
relation functions by means of differences of correlated dynamical processes 
[2, 3] (for analytic applications of these kinds of methods see for example ref. 
[4]). Here we will propose a new method which takes care about a crucial in- 
gredient, the gauge invariance of the theory, and that looks far more effective 
than the pre-existing schemes. 

The simple idea of the difference scheme to compute connected correlation 
functions is the following. In the rest of the paper we will consider the 0+^ 
glueball mass, defined from the connected correlation function 

G{t) = (^(3) (0)^(3) (t)) - {E^^\0)){E^^\t)) ~ exp{-mo++t} , (1) 

where with the upper script 3 we denote the sum of the spatial plaquettes on 
the 3-cube at a given time (0 and t in the previous equation). We will consider 
a system governed by Wilson Action at inverse temperature /? = T~^, i.e. 

Sp^P E (l-^Trf/p), (2) 

plaquettes 

for SU{2) gauge fields, with the usual definition of the plaquette variables 
Up. 

A possibly effective way to measure connected correlation functions has 
been suggested [2, 3], and it is based on simulating the dynamics of two copies 
of the system, one with the original action (2) and one with a modified action. 
In our case we only modify the Action involving links on the 3d cube ait — 
(we are in 4 euclidian space-time dimensions, and we identify one of these 
dimensions with the euclidian time). So, for t 7^ we use the original Wilson 
Action, while at t = we modify the value of the coupling, by setting 
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(3 + Sp , 



(3) 



with a small 6/3. In the limit of 5/3 one gets that [3] (we indicate with a 
tilde expectation values taken over the modified action) 



where by E"^^^ (i) we denote the average of the space-hke plaquettes at time t. 
In words we consider the difference of the time t energy of the unperturbed 
and the perturbed system, in the limit of a small perturbation. This differ- 
ence, when measured with good statistical precision, gives us a measurement 
of the connected correlation function G. 

The main point of the method is that if we use an appropriate simulation 
technique E and E arc correlated, and by exploiting that we can eliminate the 
most part of the statistical error we would get in a direct measurement of the 
G{t) . The typical pattern of such a joint simulation is the following. We start 
from two copies of the same configuration. Their energy distance is zero, since 
we are out of equilibrium. Now we start the two simulations. The energy 
distance at a given time t will hopefully stabilize after a transient period (tq, 
we indicate with r the dynamical fifth time, as opposed to euclidian time 
t) at the correct distance. Then for large dynamical time (r > ri) the two 
trajectories in phase space will separate, the statistical noise will dominate 
the signal and the measurement will be of little use. The measurement 
window goes from tq to ti, and a good dynamical procedure maximizes its 
extent. 

A discrete dynamics like the Metropolis algorithm does not do the job. 
Indeed in this case the advantage of the method, which to be well performing 
exploits the fact that the two trajectories are close in phase space, is lost due 
to the intrinsic discreteness of the updating procedure. On the contrary a 
Langevin dynamics can be the basis of our scheme [2, 3]. For the two copies 
of the system we write 



(4) 




(5) 
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where we have denoted by U the unperturbed fields and by V the perturbed 
ones, and r] is the same noise for the two systems. 

For a gauge model the gauge degrees of freedom create an additional 
complication. The gauge part of the degree of freedom random walks in 
phase space, and such a random walk tends to separate the two trajecto- 
ries. An usual gauge fixing (for example putting to one all time-like gauge 
variables) does slow down the dynamics making the method inpractical [3]. 
In ref. [3] we have seen that sometimes the phenomenon can be dramatic 
enough to make any measurement impossible (even for short euclidean time 
separations). 

In ref. [3] we proposed to solve the problem by using some kind of mag- 
netic field, which would insure a smooth, partial fixing of the gauge. Wilson 
action would be modified by a non-gauge invariant term, selected in such 
away to have a small effect one gauge-invariant quantities. Even if this 
method was improving the situation, it did not turned out to be very supe- 
rior to usual methods. Trajectories in phase space did diverge quite soon, 
and a slow drift on the system internal energy (that is also modified of a 
small amount due to the magnetic like term) was difficult to control. Also 
a non-gauge invariance formulation is definitely not so appealing, and could 
present a large number of unwanted features. 

In this note we propose a new method, and we show that it is indeed 
very effective. The basic idea is very simple. We fix the gauge where the two 
configurations are as similar as possible. We call such a gauge the chameleon 
gauge. This method does indeed keep the gauge part of the two systems 
as close as possible (reducing the rate of divergence of the two trajectories 
in phase space) but does not introduce any sizeable slowing down in the 
dynamics of observable quantities like the energy. 

After each full lattice sweep of Langevin update of the link variables U 
and V we gauge fix the U configuration (obviously it does not change to 
gauge fix instead the V configuration). We maximize the quantity 



where N runs over the lattice sites and over the 4 directions. We find the 
gauge transformation {g{n)}, with 




(6) 




(7) 
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We determine {g{n)} such that 

'^g{n)Un{n)g^n + jj,)V^{n) — local maximum . (8) 

n,n 

The quantity 

- U,{n)Vi{n)) , (9) 

links 

gives a measure of how good a gauge fix we have been able to reach. — 
when the two configurations are identical. A high value of T signals that 
a good gauge fix has not been achieved (implying for example that the two 
gauge field configuration differ because of a physical, non gauge freedom 
related reason). 

Our numerical simulations show that the method is very well performing, 
and has very pleasant features. The (correct) equilibrium value for the pla- 
quettc energy is reached in very short time. That means that all the problems 
connected to standard gauge fixing (like axial gauge) or to our magnetic field 
partial gauge fixing have been solved by the chameleon gauge. We show in 
fig. (1) the plaquette value computed with the magnetic field method, where 
a slow drift is evident. In fig. (2) we show the new data, where after less than 
100 Monte Carlo sweeps the measured value is stable. Correlation functions 
measured by difference technique on the chameleon gauged configurations 
are far less noisy than with our old method. In a typical situation we gain a 
factor of order 5 over the time window we can rely on. It is also quite inter- 
esting to note that the breakdown of the correlation functions (that become 
noisy after a given time ri) is always signaled by a sudden growth of the 
quantity i.e. by a collapse of the quality of the gauge fixing we are able 
to reach. 

Let us give a few details about our runs and our numerical results. We 
simulate an SU (2) system on a 16^ lattice with periodic boundary conditions. 
We keep one copy of the system at the original /3 value, while we simulate 
four parallel systems, where different time slices are set at the inverse square 
coupling value P + S/3. The program size is 44 Mbyte, and on a IBM RISC 
WS 550E we run 500 steps in close to one day. We have selected (3 = 2.5, 
that is close enough to the scaling region, and where our lattice size is large 
enough to describe infinite volume behavior. Michael and Teper estimate in 
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Figure 1: The plaquette operator expectation value computed by using the 
magnetic field method. 
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Figure 2: The plaquette operator expectation value computed by using the 
chameleon gauge fixing method. 
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ref. [5] that ma ~ .660. The result we obtain by using our method is fully 
compatible with that, and from a global fit to our data we obtain ma ~ .67, 
with an error, including both statistic uncertainty and systematic effects from 
finite euclidean time distance that we estimate to be smaller than 5%. 

In fig. (3) we show the signal we were getting at separation 1 by using 
our old method (based on the use of the magnetic gauge fixing term), while 
in fig. (4) we show the result we obtain with our new method. The difference 
is quite impressive. While signal was becoming noisy at r ~ 400 and was 
lost at r ~ 800 in the old approach, in the chameleon gauge fixed approach 
we get a stable signal up to time r > 1200. 

We believe that there is still much to understand in the physics of the 
so-called difference methods, both in gauge and in simple spin systems, but 
we hope that the method we are proposing here is a nice step in the direction 
of less expensive measurements of mass spectra. 
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Figure 3: The distance 1 correlation computed by using the magnetic field 
method. 



9 



Distance 1 



0.1 



0.08 - 



0.06 - 



0.04 



0.02 - 



- 



-0.02 



1 6* lattice 
P=2.5 
5p=0001 
Gauge Fix=5 



200 



400 



600 
step 



800 



1000 



1200 



Figure 4: The distance 1 correlation computed by using the chameleon 
gauge fixing method. 
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